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ABSTRACT 

We present RXTE observations of the eclipsing X-ray binary Hercules X-l conducted during 
an anomalous low state. Data reduction reveals a light curve over 2.7 orbital cycles remarkably 
similar to optical and UV light curves which are dominated by the companion star. Count rates 
are modulated close to the orbital period, attaining a maximum when the inner face of the 
companion star, irradiated by X-rays from the compact source, is most-visible. Cold reflection 
provides an acceptable fit to the energy spectrum. Employing binary geometry to scale the model 
and assuming companion star reflection, we are able to reconstruct the incident X-rays which 
are removed from our direct line-of-sight (presumably by the accretion disk). We find the flux 
of the hidden source to be identical to the observed flux of Her X-l at the peak of its main- 
high state. Consequently, Her X-l is emitting a reflected spectrum, largely uncontaminated by 
direct X-rays in the anomalous low-state. The spectral energy distribution, period, amplitude 
and phasing of the modulation are all consistent with a companion star origin. Since this source 
occurs in a well- understood binary environment, it provides an excellent case study for more 
sensitive experiments in the future. 

Subject headings: binaries: close - binaries: eclipsing - pulsars: individual (Hercules X-l) - radiative 
transfer - stars: neutron - X-rays: stars 
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1. Introduction 

Her X-l (Tananbaum et al. 1972) is an eclipsing 
X-ray binary containing a pulsar of 1.4 M Q and 
an A7 stellar companion of 2.2 M Q (Middlcditch 
& Nelson 1976; Reynolds ct al. 1997). The sys- 
tem displays behaviour on four separate periods 
- the pulsar spin (1.24-s), the binary orbit (1.7- 
d), a super period of 35-d which results from a 
retrograde-precessing, warped accretion disk, and 
a beat between the precessional and orbital pe- 
riods of 1.62-d. The warp phenomenon is poorly 
understood but its source is likely to be a combina- 
tion of both tide- and radiation-driven precession 
(Papaloizou & Tarquem 1995; Pringle 1996). 

The 35-d super-cycle results in two phases 
of strong X-ray activity per cycle (Giacconi 
et al. 1973; Scott & Leahy 1999). The main-high 
state has a rapid turn-on over ~90-m and decays 
over ^10-d. An off-state where flux remains at 
~1 percent of the main- high state level follows for 
the next ~10-d, succeeded by a second short- high 
state, lasting ~5-d, with flux peaking at ~30 per- 
cent of maximum. A further off-state completes 
the cycle and extends over the next ~10-d. The 
accretion disk casts a shadow over the surface of 
the companion star which migrates across the star 
on the beat frequency. The super-cycle is conse- 
quently observed at UV and optical energies which 
are dominated by X-ray emission reprocessed in 
the stellar atmosphere (Gerend & Boynton 1976). 

The 35-d clock has remained mostly coherent 
since it's discovery. There have however been 
three occasions when the clock has missed several 
consecutive turn-ons, or main- high flux has been 
reduced (Parmar et al. 1985; Vrtilek et al. 1994; 
Parmar et al. 1999). The cause of these anoma- 
lous low-states is not clear, but probably result 
from changes in the state of the accretion disk, ei- 
ther through an increase in vertical scale height, 
or disruption to the disk warp. However the level 
of reprocessed UV and optical flux from the com- 
panion star during these states suggests that the 
intrinsic X-ray flux, and presumably the accretion 
rate, remains mostly unaffected by this event (Vr- 
tilek et al. 1994; 2001). The high orbital inclina- 
tion of the binary indicates that perhaps just a 
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small variation in the disk geometry will block di- 
rect X-rays from our line-of-sight, while the com- 
panion stars view remains largely unchanged over 
much of its surface area. 

The latest anomalous low-state began during 
1999 Feb (Parmar et al. 1999) and the on-states 
had not recovered four cycles later when the cur- 
rent observations were scheduled. These were 
part of a simultaneous multi-wavelength campaign 
on Her X-l described by Vrtilek et al. (2001). 
Results from simultaneous Hubble Space Tele- 
scope orbits are detailed by Boroson ct al. (2000) 
and Boroson et al. (2001), and an earlier epoch 
of short-high state visits with RXTE by Still 
et al. (2001). 

The current RXTE data reveal many similar- 
ities between the X-ray light curve and its opti- 
cal/UV counterparts, suggesting that much of the 
anomalous low state X-ray emission could be of 
companion star origin. We find that cold Comp- 
ton reflection is a statistically acceptable model 
for the energy spectrum and a reflection surface 
the size and distance from the pulsar of the com- 
panion star gives a good prediction of the incident 
X-ray flux provided this is identical to the peak of 
the normal on-state intensity. The possible detec- 
tion of a relatively clean reflection spectrum in a 
well-determined binary geometry, free of dominant 
direct X-rays, would be a valuable resource for 
investigators of reflecting atmospheres in the less 
well-determined environments of AGN and stellar 
black hole candidates (Lightman & White; White, 
Lightman and Zdziarski 1988). 

2. Observations 

RXTE pointed at Her X-l intermittently be- 
tween 1999 July 11-16 (MJD 51370.83-51375.37) 
accumulating 100-ksec of time-tagged events. 
Event reconstruction was performed using stan- 
dard algorithms within FTOOLS v4.2. We analyze 
data from the Proportional Counter Array (PCA; 
Zhang ct al. 1993) in which three of the five identi- 
cal Xe Proportional Counting Units (PCUs) were 
operational with a combined effective area of 3900 
cm 2 . As well as the standard data formats, we em- 
ploy two event analyzers in the GoodXenon event 
mode with 2-s readout - timestamps are resolved 
to l-/is, while employing the full 256 energy chan- 
nels (2.2-119.0 keV during epoch 4). Response 
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Fig. 1. — The 3-30 keV band light curve over the duration of the observations. The vertical dashed 
lines correspond to superior conjunction of the neutron star, according to the orbital ephemeris of Deeter 
ct al. (1991). 



matrices from 1999 Mar 30 were obtained from 
the HE AS ARC archive 3 . Background estimates 
are derived from the Very Large Event model 4 to 
account for cosmic events, internal particle gen- 
eration and South Atlantic Anomaly activation. 
Events from the High Energy X-ray Timing Ex- 
periment (HEXTE), also onboard the spacecraft, 
proved too scarce to be useful. 

Throughout this paper we will employ the or- 
bital ephemeris of Deeter et al. (1991): 

T orb = MJD 43 804.51998(1)+ 1.70016772(l)£ orb 

-5.2(5) xlO- 11 ^! (1) 

where T orb corresponds to superior conjunction of 
the compact star and E orb is the orbital cycle num- 
ber. Assuming the 35-d epoch is determined by: 

T 35 = MJD 50 041.0 + 34.85 E 35 (2) 

where T 35 corresponds to X-ray turn-on of the 
main-high state and E 35 is the cycle number 

3 http:/ /xtc. gsfc.nasa.gov 

4 http://lheawww. gsfc.nasa.gov/~stark/pca/pcabackest. html 



(M. Kunz, private communication), the precession 
phases sampled are </> 35 = 0.16-0.29. These would 
normally correspond approximately to the the sec- 
ond half of a main-high state (Scott and Leahy 
1999), but the average count rate of 4 counts 
s _1 PCU -1 testifies to the anomalous extended 
low-state in which Her X-l had remained for the 
previous four 35-d cycles (Coburn et al. 2000). 
Count rates are of order 1 percent of the normal 
main-high state (dal Fiume et al. 1998). 

3. Light curve 

Standard2 data from each visit were filtered 
to reject pointings closer than 10 degrees to the 
Earth's limb and those off-axis by 0.02 degrees or 
greater. From all PCA layers and columns, events 
were summed across pulse-height channels to pro- 
vide a light curve in the energy range 3-30 keV 
with 1-ksec sampling. This is presented in Fig 1 
and hardness ratios between 2-5 keV, 5-8 keV and 
8-20 keV are provided in Fig 2 (eclipse data have 
been removed from the hardness diagrams). Back- 
ground models have been subtracted and data di- 
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Fig. 2. — Hardness ratios for the bands 2-5 keV, 5-8 keV and 8-20 keV. Eclipse points, <fi orb 
have been excluded. 
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vided by the number of active PCUs. 

Gerend & Boynton (1976) show that the major- 
ity of power from companion star emission occurs 
on the 1.7-d orbital frequency through stellar ro- 
tation. However, light curve structure associated 
with the accretion flow from the companion star, 
such as the hot spot and splash from the colli- 
sion of stream and disk (Schandl 1996), occur on 
a timescale close to the beat period. Fig 1 shows 
an intensity modulation approximating the 1.7-d 
orbital period. Depending on its source, we ex- 
pect the observed modulation to occur cither on 
the orbital or 1.62-d beat period. 

The hardness ratios of Fig. 2 indicate that the 
modulation is not the result of photo-electric ab- 
sorption. We are observing intrinsic flux increases 
or opaque obstruction of a soft emission region. 
The observed structure occurs over timescales sev- 
eral factors longer than normal main-high dipping 
events (Leahy 1997). 

Maximum flux corresponds approximately with 
orbital phase (f) OIb = 0.5, where </> orb = 0.0 is su- 
perior conjunction of the neutron star. This is 



consistent with a region associated with the irra- 
diated inner surface of the companion star. The 
similarities between this light curve and the op- 
tical, UV and EUV light curves which are dom- 
inated by companion star emission are striking 
(Still et al. 1997; Leahy & Marshall 1999; Boro- 
son et al. 2000). We suggest that the modulation 
is the result of Compton reflection of X-rays off 
the atmosphere of the companion, which was pre- 
dicted by Basko, Sunyaev & Titarchuk (1974) and 
detected at EUV wavelengths by Leahy & Mar- 
shall (1999) and Leahy et al. (2000). 

If we are detecting companion star reflection, 
there must be a second source of X-ray flux asso- 
ciated with the compact object in order to produce 
the narrow eclipses at MJD 51 372.0, MJD 51 373.7 
and MJD 51 375.4. There are equivalent eclipses at 
optical and UV wavelengths which are thought to 
be the eclipse of disk light by the companion star 
(Gerend & Boynton 1976; Boroson et al. 2000). 

After removing eclipses and using the Lomb- 
Scargle statistic (Scargle 1982), the best fit period 
to this light curve is 1.74 ± 0.06-d, with a false 
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Fig. 3. — A fit to the light curve using the geometric binary synthesis code of Still et al. (1997) plus a 
linear-decay function to represent the corona. Points designated by open squares, corresponding to eclipses 
and the dip at MJD 51,371.05 were ignored by the fit. The dashed line is the coronal fit alone. The three 
spectra presented in Fig. 4 were extracted from the visits indicated by arrows. 



alarm-probability of 10~ 12 . The error is the a- 
width of the peak in the power spectrum and can 
only discriminate between the orbital and beat fre- 
quencies to 1.3-cr. 

We attempt a crude fit to the light curve where 
the reflection component is constructed by com- 
bining the geometric stellar synthesis code of Still 
et al. (1997) with the classical H-function reflec- 
tion coefficients (Chandrasckhar 1960). Reflected 
flux between 3 and 30 keV is integrated over 
the stellar surface area. This is detailed in Ap- 
pendix A. We adopt an orbital inclination of i 
= 82°, a neutron star mass of M ns = 1.4 M Q 
and a Roche lobe-filling companion star mass of 
M com = 2.2 M Q (Reynolds et al. 1997). Unlike 
Still et al. (1997), and in order to reduce the num- 
ber of free parameters, we do not attempt to model 
the shadow of the accretion disk over the compan- 
ion star, or the orbital eclipse of the star by the 
disk, centered at </> orb = 0.5. The energy spectrum 
of the compact source is taken to be a powerlaw 
of slope a = 0.9 with an exponential cutoff at 
E cut = 11 keV (see Sec. 4). The luminosity of the 



source is arbitrary, where the resulting reflection 
spectrum is rescaled such that the maximum count 
rate at orb = 0.5 is 1 count s _1 PCU -1 . The coro- 
nal component is represented by a linear function 
which is better statistically than a constant, such 
that the overall fit is given by: 

L(t) = Pl + p 2 t + P3 R(t) (3) 

Pi and pi are the constant and linear terms of 
the coronal component, and p 3 is a renormaliza- 
tion constant for the modeled companion star light 
curve R(t). t is time in days after MJD 51,000. 

After rejecting eclipse points and the dip at 
MJD 51 371.05, the fit was made using a Downhill 
Simplex Algorithm (Press et al. 1992). The result 
is plotted in Fig. 3 where p\ = 145, p 2 = —0.38 
and p 3 = 2.14; \ 2 = 2100 for 84 dof. Decreas- 
ing coronal flux is consistent with the behaviour 
of direct photons from the accretion region during 
this phase of the normal 35-d cycle (Scott & Leahy 
1999). 

The fit is poor, although this is not particu- 
larly discouraging given our first-order treatment 
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of coronal variability. Residuals are probably dom- 
inated by the coronal component. In the follow- 
ing sections we perform spectral fits in an attempt 
to separate the two components, test the suitabil- 
ity of reflection for fitting the energy spectrum of 
the modulated component, and subtract modeled 
stellar reflection spectra to reveal the spectral vari- 
ability of the corona. 

4. Energy spectrum 

Under the same filtering constraints as above, 
the GoodXenon events were sampled with the full 
pulse-height resolution of the PCA. Channels 1- 
7 (< 3 kcV) were ignored due to uncertain back- 
ground modeling, and channels 63-256 (> 30 keV) 
ignored as a result of poor counting statistics. 
xspec vll was used at all times. Spectral fits 
to three of the visits are presented in Fig 4. 

4.1. Corona spectrum 

We consider a model where emission originates 
from two distinct binary regions, an accretion disk 
corona and the companion star. We first consider 
the corona, whose output is presumably the result 
of upscattering soft disk photons in a relatively 
hot medium (e.g. Sunyaev and Titarchuk 1980). 
For now, the parameters of this model will be con- 
sidered constant over time (excluding eclipses), ex- 
cept for the normalization parameter which we will 
vary according to the linear component of our light 
curve fit from Sec. 3. This rather false constraint 
will allow us to separate approximately the two 
spectral components and is probably valid across 
neighboring visits, but not over a number of or- 
bital cycles (Fig. 3). 

The visit at MJD 51372.15 is the most likely 
to isolate the coronal component. Eclipse egress 
is approximately complete and any contribution 
from the reflection component is assumed negligi- 
ble. We refrain from using the pre-eclipse visits at 
MJD 51 373.4 and MJD 51 375.2 in case dipping is 
occurring. An absorbed power law model of the 
form: 

S(E) = cxp [-N H a(E)}AE- a (4) 

where A is the normalization, <r(E) the photo- 
electric cross-section (Balucinska-Church and Mc- 
Cammon 1992) and N H the neutral hydrogen col- 
umn density, provides an inadequate fit of x 2 od = 



2.06 for 54 dof. The addition of a Gaussian model: 
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representing an iron emission complex at 6.4-6.7 
keV improves the fit to x 2 od = 1-15 for 51 dof. 
A K is the line flux, cr K the line width and E K the 
line center. Compared to previous observations 
of Her X-l in anomalous low-states (e.g. Par- 
mar et al. 1999; Coburn et al. 2000), the pow- 
erlaw component is rather steep, a — 1.8 ± 0.1, 
and the Fe line nested at high energy 6.9 ± 0.2 
keV. By introducing further complexity into the 
model both the powerlaw slope and line energy 
agree with previous results but the x 2_s tatistic 
cannot discriminate between a number of models. 
A definitive model is not essential since these fit 
parameters will remain fixed in order to extract 
approximately the reflection spectrum a few visits 
later. Therefore we choose a model which is sim- 
ple but physically realistic, and is also compatible 
with fits made earlier in the anomalous low state 
by Coburn et al. (2000) - in this case a powerlaw 
with high-energy exponential cutoff, a Fe emission 
line and a partial-covering, cold absorber: 



S(E) = P{E) 
where 



.loxp ( ~^-)E " rG[E) 



(6) 



P(E) = (l-f) + fexp[-N H a(E)) (7) 

/ is the partial covering fraction and E cat the cut- 
off energy. The excessive number of free parame- 
ters mean that the column density, covering frac- 
tion and powerlaw slope are poorly constrained. 
Since the current anomalous low-state is thought 
to be the result of obscuration rather than intrin- 
sic source behaviour we can assume that the coro- 
nal photons have undergone elastic scattering and 
rationalize fixing the powerlaw slope to it's main- 
high state value, a = 0.9 (Dal Fiume et al. 1998). 
From the visit at MJD 51 372.15, N H = (1.6 ±0.8) 
xlO 23 cm" 2 , / = 0.7 ± 0.1, E CM = 11 ± 2 keV, 
A = (2.6 ± 0.9)xl0~ 3 photons cm" 2 s" 1 keV^ 1 
at 1 keV, E K = 6.7 ± 0.1 keV and a K = 0.7 ± 0.2 
keV with a line strength of (2 ± l)x 10 -4 photons 
0.85 for 50 dof and the equiva- 
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lent width of the Fe complex is 775 eV. 
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Fig. 4. — Spectral fits from three visits to Her X-l. The top left panel is a powerlaw + Gaussian profile fit to 
mid-eclipse data (MJD 51371.98). The right hand panel presents the best fit model. The middle two panels 
are the corresponding fit and model of the corona for the post-eclipse visit at MJD 51372.15 consisting of 
an exponentially-cutoff powerlaw and line, absorbed by a partial cover. The lower panels represent the visit 
at maximum flux, MJD 51372.90, where a reflection component and second line are added. See Sec. 4 for 
details of the spectral models. 
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We apply this same model to a visit at (f> oth = 
0.5 (MJD 51 372.90). Keeping a and £ cut fixed at 
their previous values provides a poor fit, xiL = 
2.74 for 51 dof. Relaxing both these constraints 
yields an acceptable fit of x 2 od = 0.74 for 49 
dof. However N H = and a is extremely flat, 
0.13 ± 0.12. Consequently, in a time-variable ca- 
pacity the coronal model alone does not fit the 
data sensibly. Therefore a two-source model is 
more appropriate. 

4.2. Reflection spectrum 

We would like to test whether reflection pro- 
vides a good fit to the modulated component. Re- 
flection would be the result of scattering of hard 
radiation from the accreting source off the rela- 
tively cold atmospheres of the companion star or 
accretion disk. The spectrum will suffer photo- 
electric absorption at energies less than 10 keV, 
have a strong Fe Ka line and will be down- 
scattered at energies higher than 50 keV. Thus, the 
reflection spectrum has a maximum in the RXTE 
PCA band at energies 15 - 20 keV. 

As explained in the previous section, we assume 
the coronal model parameters are time-invariant, 
except for the normalization which is scaled by the 
light curve fit of Eqn. 3. This is in order to reduce 
the number of free parameters which is excessive 
for the quality of spectra from individual visits. 

The visit we consider is at (j) orb ~ 0.5 (MJD 
51 372.90) which has a number of advantages over 
other observations. This phase provides the best 
counting statistics. Furthermore we expect to be 
observing <~ 100 percent of the irradiated surface 
at this time so we can assume that all reflected 
flux is recorded. Finally we are reasonably close 
to the visit where the coronal spectrum was char- 
acterized and so our assumptions concerning the 
form and brightness of this component are hope- 
fully valid. We test these assumptions in Sec. 6. 

We first attempt adding a few simple models to 
the coronal component in an attempt to produce 
a statistically-acceptable fit. A second Gaussian 
component is always required. Adopting S(E) + 
power law + Gaussian provides a x 2 od = 1.17 for 52 
dof. The photon index is flat, a — 0.33 ± 0.1, and 
the fit yields a notable low-energy excess. A S(E) 
+ blackbody + Gaussian yields x 2 ed = 0.72 for 52 
dof. The blackbody temperature is 6.2 ± 0.7 keV 



and the 6.4 ± 0.1 keV line has an equivalent width 
of 3.13 keV. Despite the good fit, it is not ob- 
vious what would radiate thermally at this tem- 
perature within the framework of the anomalous 
low-state model (Vrtilek et al. 2001), or produce 
such a broad line. The low-energy excess remains 
a feature. 

More realistically, we fit S(E) + Gaussian + the 
clastic, cold Compton reflection model hrefl. Al- 
though the geometry of this model is more suited 
to disk reflection than a Roche surface, it is the 
most-appropriate model within the xspec pack- 
age and yields an approximate fit. This provides 
confidence in a more sophisticated approach un- 
dertaken in Sec. 6. 

The fraction, / esc , of direct X-rays observed is 
zero, since the point source is obscured by the 
accretion disk. However much of the companion 
star atmosphere has a largely-unobscured view of 
the X-ray source. The reflection model assumes 
an isotropically-radiating point source of a given 
energy spectrum and flux that irradiates a struc- 
ture subtending a solid angle to on the sky with 
respect to the point source. By adopting stellar 
masses of M ns = 1.4 M and M com = 2.2 M Q , and 
companion radial velocity of 100 km s _1 (Reynolds 
ct al. 1997), we find u> — 0.15-7T. We approximate 
a reflecting structure incident to X-rays at all an- 
gles between 6* min = degrees and # max = 90 deg, 
weighted equally, and inclined on the sky at an an- 
gle 9 ob3 = deg to the observer at (p otb = 0.5 and 
orbital inclination i ~ 90 degrees. Solar abun- 
dances were adopted. The incident energy spec- 
trum is modeled by the same cut-off powerlaw as 
before with a = 0.9 and E cut = 11 keV. 

The photon flux of the irradiating powerlaw 
spectrum and Gaussian profile provide just four 
free parameters, yielding a x 2 od of 0.89 for 53 dof. 
Fit parameters are listed in Tab. 1. The fit is 
extremely good considering there is only one free 
parameter available to define the reflection contin- 
uum. Encouragingly, photon flux for the unseen 
powerlaw source at 1 keV is (9.0±0.9)xl0~ 2 pho- 
tons cm~ 2 s _1 keV -1 . The flux from a typical 
main-high state of Her X-l is 0.1 photons cm~ 2 
s" 1 keV" 1 at 1 keV (dal Fiume et al. 1998). The 
line is strong, narrow and nested at low Fe Ka en- 
ergies, as predicted by reflection models (Basko 
et al. 1974; Ross & Fabian 1993). This consis- 
tency, combined with the period and phasing of 
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the modulated component all point to the com- 
panion star providing the reflective surface. 

4.3. Eclipse spectrum 

The visit during mid-eclipse (MJD 51 371.98) is 
fit well with a single powerlaw, x 2 od = 0.93 for 
26 dof over channels 8-35 (3-15 keV). Although 
the inclusion of a Fe line is not formally required, 
there is an excess longward of the Fe K edge. In- 
cluding a Gaussian model suggests a line of equiv- 
alent width 2.6 keV, where xL = 0.39 for 23 dof. 
Here a = 1.2 ± 0.3, A = (2 ± l)xl0~ 4 photons 
keV- 1 cm" 2 s- 1 at 1 keV, E K = 6.7 ± 0.3 keV 
and cr K = 0.3 ± 0.2 keV with a line intensity of 
(6±3)xl0~ 5 photons cm -2 s _1 . Adopting a cut- 
off powerlaw model of a = 0.9 and E cat = 11 keV 
and fitting just the normalization and line pa- 
rameters yields consistent Gaussian characteristics 
within the uncertainties. 

5. Pulse timing 

In this section we report the detection of spin 
pulses, measure their frequency and present pulse 
profiles. GoodXenon data were extracted over 
all energy channels and binned to 0.02-s. Times 
were corrected to the compact sources local stan- 
dard of rest using the orbital time-delay semi- 
amplitude determined by Deeter et al. (1991) and 
the spin period of the pulsar was searched for 
using the Lomb-Scargle algorithm. A period of 
1.237746±0.000001-s was found with a false-alarm 
probability of 10~ 5 . This is consistent with the 
measurements of Parmar et al. (1999) and Coburn 
et al. (2000) during the same anomalous low-state 
and signifies that the pulsar has spun down dur- 
ing the low-state episode, as it has in the previ- 
ous two recorded anomalous low epochs (Vrtilck 
et al. 1994). Coburn et al. (2000) discuss the im- 
plications of these spin-down episodes with regard 
to the accretion torque in the system. The pulse 
profile is presented in Fig. 5. 

It is unlikely that the observed pulses are within 
reflected light from the companion star. The spin 
frequency is constant in the rest-frame of the com- 
pact object so should be scattering off the accre- 
tion disk or accretion disk corona. Pulse fractions 
are of the order a few percent whereas models of 
reflected pulses off the companion predict fractions 
an order of magnitude less (Middleditch and Nel- 



son 1976). 

6. Modeling Roche lobe reflection 

The current family of reflection models avail- 
able in XSPEC Vll are not well-suited to our prob- 
lem. These were developed for accretion disk re- 
flection problems (see e.g. Magdziarz & Zdziarski 
1995) and limit our results because of incorrect 
geometry. The models available assume the re- 
flecting body is a uniform slab, viewed at an in- 
clination angle and illuminated by a point source 
above the surface of the slab. The size of the slab 
and the distance from the point source are char- 
acterized by a filling factor - the fraction of sky 
occupied by the slab with respect to the source. 
The model adopted in Sec. 4.2 takes into account 
a range of equally-weighted irradiance angles to 
approximate the curved Roche surface of the star, 
but it is again a simplification of the problem (al- 
though a valuable consistency check for more so- 
phisticated models). The simplified model did in- 
dicate however that the reflected flux was identical 
(within the 1-a measurement uncertainty) to that 
expected from a point source with the distance, 
energy spectrum, and intensity of Her X-l at the 
peak of its main-high state. We would like to test 
this more rigorously by defining an accurate re- 
flection surface. 

In order to treat the problem more accurately 
we develop a code which divides the surface of 
a Roche lobe-filling star into small, discrete ele- 
ments. Each has an individual inclination angle 
with respect to the Earth. Each has an individ- 
ual grazing incidence angle and filling factor with 
respect to the X-ray source. Our task is to de- 
termine these parameters over a grid using Roche 
geometry and sum the reflection spectrum over all 
observer- visible surface elements. This approach 
both treats the geometry correctly and allows a 
fit to be made at any orbital phase. Again, we 
follow the algorithm described in Appendix A. 

The model was developed as an external, mul- 
tiplicative XSPEC package and applied to the spec- 
trum sampled at MJD 51372.90 where we fit the 
composite model of Eqn. 6, by replacing the hrefl 
component and adopting the fixed coronal and 
pulsar parameters listed in Tab. 1. Our reflection 
parameters are stellar mass ratio q = M com /M na , 
which defines the shape of the companion star 
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Table 1 

best-fit parameters for the three spectra presented in flg. 4. all fits are also subject 
to a fixed interstellar column density of 5.1xl0 19 cm~ 2 (dal flume 1998). the reflection 

parameters are given for mjd 51 372.90. 



Eclipse parameters 



a 1.2 ±0.4 

A (ph cm- 2 s- 1 kcV- 1 ) (2 ± 1) x 10~ 4 

E K (kcV) 6.7 ±0.3 

cr K (kcV) 0.3 ± 0.2 

A K (ph cm~ 2 s' 1 ) (6±3)xl0~ 5 

Coronal parameters 

iV H (xHP cm" 2 ) 16 ± 12 



/ 0.7 ±0.6 

a a 0.9 

E cM (keV) 11 ± 2 

A (ph cm" 2 s- 1 keV" 1 ) (3 ± l)xl0~ 3 

E K (kcV) 6.7 ±0.2 

<T K (keV) 0.7 ±0.4 

A K (ph cm~ 2 s' 1 ) (2± I)xl0~ 4 

Reflection parameters 

a a > 6 09 

E cut (keV) a <" 11 

A (ph cm" 2 s- 1 keV- 1 ) 13 (9.0 ± 0.9) x 10~ 2 

U & 0.157T 

cos 6> oba a 1 

cos 6> mi „ a 1 

cos 6> max a 

Fe abundance (solar) a 1 

f a 

J esc w 

£ K (keV) 6.5 ±0.1 

cr K (kcV) 0.1 ±0.1 

A K (ph cm" 2 s- 1 ) (1.5±0.2)xl0~ 4 

a Frozen parameter. 

b Model component is directly-invisible to the 
observer. 
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Fig. 5. — Time-averaged spin pulse from Her X-l. The pulse cycle is repeated for clarity. 



and the area of pulsar-sky each surface element 
occupies, the orbital inclination, i, and the or- 
bital phase, <p otb , which together define those el- 
ements of the stellar surface visible to Earth. We 
use the spectroscopic measurements of Reynolds 
et al. (1997), q = 1.57 and i = 82°. The orbital 
ephemeris of Eqn. 1 indicates 4> orb — 0.55. In or- 
der to calculate the scattering probability the Fe 
abundance is adopted as solar. 

As before, four free parameters remain, yield- 
ing a reduced % 2 of 0.86 for 53 dof, where A = 
(9.6 ± 0.5) xlO -2 photons s" 1 cm~ 2 keV" 1 at 1 
keV, E K = 6.5 ± 0.1 kcV, a K = 0.1 ± 0.1 keV 
and A K = (1.4 ± 0.2) xlO" 4 photons s" 1 cm" 2 . 
This is both consistent with the hrefl model and 
confirms that the intrinsic flux from the pulsar is 
reproduced with reflection off a Roche surface at 
this orbital phase. 

At <j> OIb = 0.5 the Roche reflection model most- 
resembles the disk reflection models. Ideally we 
would like to apply the test above to all orbital 
phases but Fig. 3 indicates that if stellar reflec- 
tion drives the orbital light curve, then the coronal 
spectrum is variable and our model assumptions 
break down. The alternative approach is to freeze 
A, calculate the reflection spectrum at each visit 
and then fit the coronal model to the residuals. 

Since it is not possible to separate the coronal 
and reflected Fe lines, these are fit with just a sin- 
gle Gaussian. Neglecting the short-term behaviour 
of the absorbing column and coronal line flux, we 
expect Gaussian strength, Ak, and width, ok, to 
vary with the reflection continuum as a result of 
the narrow fluorescence line from the companion 



star varying over the orbital cycle. The 3-30 keV 
reflection flux (top panel) and the best fit parame- 
ters for all visits are presented in Fig. 6. Both the 
coronal a and E cM were frozen during the fits. 

Coronal parameters are expected to vary on the 
35-d, rather than the 1.7-d cycle and we see no 
detectable orbital coherency in N H , f or A. The 
decay of the coronal flux appears related to a long- 
term increase in N H , rather than a decrease in A, 
which is consistent with the 35-d disk precession 
model for Her X-l (c.f. Scott & Leahy 1999). As 
predicted, Ak is modulated on the orbital period, 
with a maximum at <p orb ~ 0.5, although similar 
coherency is not discernible in uk- 

7. Discussion 

Despite the reflection effect in Her X-l being 
observed in the EUV band (Leahy & Marshall 
1999; Leahy et al. 2000), this is the first convinc- 
ing detection of the phenomenon at X-ray ener- 
gies. Basko et al. (1974) demonstrated that X-ray 
reflected flux from the companion star in Her X-l 
should be at the level a few percent of the intrin- 
sic radiation. This should generally be washed out 
by shot noise from the central source during nor- 
mal main-high states, but as in the current obser- 
vations, should also be observable during normal 
low states of the 35-d cycle. Sheffer et al. (1992) 
claimed a tentative detection of reflection during a 
normal off-state during the ASTRON mission, but 
a systematic search for this effect during normal or 
anomalous low-states has never been undertaken; 
the HEASARC archive is not rich with data from 
these epochs and those that exist coincide with 
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eclipses of the pulsar (e.g. Choi et al. 1997). Ex- 
ceptions are the EXOSAT pointings by Parmar 
et al. (1985), the BeppoSAX observations of Par- 
mar ct al. (1999) and the RXTE visits by Coburn 
et al. (2000). None of these analyses consider re- 
flection as a possible interpretation of data, but 
the count rates presented by Parmar et al. (1985; 
1999) arguably both display 1.7-d modulation. 
The baseline of the Coburn et al. (2000) visits is 
20 ksec covering <fi orb = 0.64-0.77. The general 
trend in count rate over this period is one of de- 
cline, as predicted by the reflection model (Fig. 3), 
although the count rate and decline are a factor 
two larger than the current data. 

Potentially, Her X-l contains a rich lode of re- 
flection physics. The reflection effect is common 
to a large number of different X-ray sources such 
as AGN and stellar mass black hole (Lightman & 
White; White et al. 1988). But generally there is 
uncertainty concerning the physical location and 
size of hard radiation sources in these objects, and 
the size and shape of the reflecting atmosphere 
(which in these cases is presumably an accretion 
disk). Therefore Her X-l is an ideal X-ray source 
for studying the properties of reflection and asso- 
ciated processes (such as winds and coronae etc) 
because we have a relatively clear picture of com- 
panion star geometry and its illumination by well- 
measured hard radiation from a central object. 

Despite this there are some limitations. Firstly 
it is not clear how, or whether, the X-ray emission 
from the pulsar is beamed, or how column den- 
sity changes with disk inclination angle. Although 
the high inclination of the binary means the ob- 
server has a similar line of sight to the source as the 
equatorial regions of the companion star, surface 
elements closer to the stellar poles may perceive 
the central source differently. In this paper we 
have assumed the source is an isotropic emitter. 
Also we assume that the X-ray source is point- 
like. This is not true, since we observe scattered 
light from a disk corona and expect further scat- 
tering and reflection from the disk itself; Leahy 
(1995) showed that soft photons come from an ex- 
tended region using eclipse timings. Furthermore 
we made no attempt to model the accretion disk 
shadow, or its eclipse, over the equatorial regions 
of the companion star since the disk shape is so 
poorly constrained. The reflection process itself is 
also treated simply. More sophisticated models by 



(e.g. Ross & Fabian 1993; Nayakshin 2000) con- 
sider the physical response of an atmosphere to 
irradiation and the resultant spectral changes. 

There is scope for much more observation. The 
reflection model can be tested more sensitively 
with improved models and greater energy resolu- 
tion using grating spectrometers sensitive to soft 
energies below 1 keV which should be populated 
by resonance absorption edges (Baskoet al. 1974), 
and over the Fe a complex to resolve and separate 
line components between corona and companion 
star. In order to constrain the location of reflec- 
tion, timing experiments with improved sampling 
over the 35-d cycle would determine whether light 
curve modulation is driven on the orbital (1.7-d) 
or beat (1.6-d) period. Fast timing visits during 
eclipses of the companion by the accretion disk at 
0orb = 0.5 can help constrain the shape of the ac- 
cretion disk using the eclipse ingress and egress 
profiles. 

8. Conclusion 

Over 2.7 orbital cycles during the anomalous 
low-state of 1999 the X-ray pulsar Her X-l dis- 
played an X-ray light curve very similar to optical 
and UV counterparts, suggesting that the com- 
panion star contributes to X-ray emission. A cold 
Compton reflection spectrum was successfully fit 
to the variable component of the energy spectrum. 
By employing binary geometry and assuming com- 
panion star reflection, the intensity of the hidden 
pulsar was found to be identical to the main-high 
state of the source when the pulsar is visible. This 
result indirectly implies the anomalous low-state 
model is correct - intrinsic accretion luminosity 
remains constant, but is occulted along the ob- 
servers line-of-sight by the accretion disk. 
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01-ADP-108, NAG5-6711 and NAG5-7333. KO, 
KH and HQ acknowledge support from the Par- 
ticle Physics and Astronomy Research Council. 
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the High Energy Astrophysics Science Archive Re- 
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Goddard Space Flight Center. 
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A. Classical Compton reflection off a Roche lobe-filling surface 

In order to estimate the stellar-reflected flux one must know the fraction of X-rays intercepted by the 
atmosphere of the companion and also the reflection properties of that atmosphere. The following calculations 
construct a stellar surface confined by the Roche cquipotential surface of a binary of stellar mass ratio 
q = M com /M ns (Sec. A.l). They then determine the monochromatic Compton reflected fraction off discrete 
surface elements and sum the reflected flux in those elements visible to the observer at binary inclination i 
and binary phase </> orb (Sec. A. 2). 

A.l. Roche geometry 

The critical Roche potential is defined by: 

* cr= <L + i + (l±^t (A1) 

r cr 1 — r cr 2 

r cr is normalized to the stellar separation and determined by minimizing \e CT /(de CI /dr CI )\, where 

e cr = (l + g )r cr -l- A + _L_ (A2 ) 

Using a spherical polar coordinate system {r,9,4>), with an origin at the center of mass of the companion 
star, the stellar surface is divided over 9 and 4> into segments of equal solid angle (sec Fig. 7). For any given 
9 and <f), r is determined by minimizing |e/(de/dr)|, where 

e = 2 + (1 + r 2 - 2rcos9) 1 / 2 + (1 + ^ C ° S - + sin 2 9 cos 2 </> + - 1 - - rcostf - * cr (A3) 
r 2 2(1 + q) 

Following the approach previously adopted by Rutten and Dhillon (1994) and Still et al. (1997), each 
surface segment is divided into two, such that the critical Roche surface is completely tiled by N triangles. 
This provides a better approximation to a curved surface than other polygons. The coordinates of each 
triangle vertex, ~Vj k = ( r jk,Sjk,4>jk), where j = 1 to 3 and k = 1 to N, are calculated using Eqns. A1-A3 
and converted to a Cartesian coordinate system (xj k ,yj k , z jk) with an origin at the center of mass using: 

Xjk = r jk cos 9 jk - z ~ (A4) 
1 + q 

Vjk = rjk sin 9 jk cos 4> jk (A5) 

Zjk = rjk sin 9 jk sin <f> jk (A6) 

The coordinate center of each triangle is = J2^j = i( x jk/3,yj k /3, Zj k /3). The unit vector normal of each 
triangular surface is 

(V 2k - Vi fc ) x (V 3k - V lk ) 

n k — — — (A7) 

\(V 2k - V lk ) x (V 3k - V lfc ) | V ; 

where the surface area of each element A k = \n k \ /2. The cosine of a surface elements normal vector and 
the unit direction vector to Earth, e, is the incidence angle of reflected photons: 

A*fc = n k .e (A8) 

e is defined using the orbital inclination of the binary, i, and the orbital phase, <p OTb . 

e x = sin i cos - 2<p orb )) (A9) 

e y = sinisin(7r(l - 2</> orb )) (A10) 
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e z = cos i 



(All) 



Similarly the cosine of each elements normal vector with the unit vector direction to the X-ray source yields 
the incidence angle of photons arriving at the stellar surface: 



n k . 



\dh\ 



where 



<7 



1 + 9 



(A12) 



(A13) 

(A14) 
(A15) 



and \dk\ is the distance between the source and each triangle. We assume the X-ray source is point-like and 
is coincident with the center of mass of the pulsar. Finally the solid angle of sky subtended by each triangle 
with respect to the X-ray source is 



cos tan" 



y/A k /ir 
\d k \ 



(A16) 



A. 2. Reflected flux 



In Chakrabarti and Titarchuk (1995) the reflection problem in a plane geometry was solved exactly using 
the Fokker-Plank treatment for multiple scattered photons. Photoelectric absorption and downscattering 
effects were taken into account correctly. They suggest that in the very non-relativistic energy range, 
< 30 keV the downscattering effects are negligible. Photons scatter off electrons almost coherently and 
consequently the results of the classical reflection problem are applicable (Chandrasekhar 1960). If the 
intrinsic flux spectrum from a point source radiating isotropically from the location of the pulsar is Fe, the 
reflected flux observed off all surface elements is 



!e = V ^^H(\ E , LL k )H(\ E , Vk ) (jM k >0,v k > 0) 

k~i^ k Vk 



(A17) 



Unobservable elements with [i k < and elements not incident to photons from the pulsar, v k < 0, should be 
removed from the sum. A# is the photon scattering probability approximated by the piecewise polynomial fit 
of Morrison & McCammon (1983). H(Xe, (3) is the classical H-function (Chandrasekhar 1960), approximated 
with an uncertainty of < 1% using the analytical expression of Basko (1978): 



H(X E ,I3) 



1 + V3/3 



1 + v/3(l - X E )P 



XeP 



(1 + A|)(ln/3 + 1.33 - 1.458/3 62 ) 



(A18) 
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Fig. 7. — Schematic of the geometry used to determine the reflected spectrum, both cartesian and polar 
coordinates are displayed, v is the incidence angle of reflection relative to the neutron star. The incidence 
angle between surface elements and the observers line of sight, \i, is the projection angle of each surface 
element on the page. 
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